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Abstract 

We briefly discuss the recent developments in causal dissipative hydrodynamic for relativistic 
heavy ion collisions. Phenomenological estimate of QGP viscosity over entropy ratio from several 
experimental data, e.g. STAR's <j) meson data, centrality dependence of elliptic flow, universal 
scaling elliptic flow etc. are discussed. QGP viscosity, extracted from hydrodynamical model 
analysis can have very large systematic uncertainty due to uncertain initial conditions. 
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1. Introduction 

In recent years there is much interest in QGP viscosity. In any fluid motion, some 
amount of irreversibility is present that leads to energy dissipation. It is then important 
to have estimates of dissipative coefficients like viscosity of QGP. Without this estimate, 
it is difficult to claim that the matter produced in RHIC collisions can be understood, 
at least in part, with thermodynamic concepts. Theoretically, it is possible to estimate 
dissipative coefficients of QGP from first principle, however, the problem is complex and 
yet unsolved. One hopes to extract viscosity from experimental data. It require numer- 
ical implementation of dissipative hydrodynamics. In recent years significant progress 
has been made in numerical implementation of dissipative hydrodynamics [1-9]. In the 
following, we briefly describe causal dissipative hydrodynamics. We also discuss about 
some estimates of shear viscosity to entropy ratio of QGP from experimental data. 
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2. Causal dissipative hydrodynamics 



Relativistic generalisation of Navier-Stokes equation, called 1st order theories [10,11] 
for dissipative hydrodynamics, suffer from the problem of acausality and instability [12]. 
The problem is eliminated in Isracl-Stcwart's [13] 2nd order theory (for detailed exposi- 
tion see [14,15]). Briefly, dissipative fluxes are treated as extended thermodynamic vari- 
ables and relaxation equations for the dissipative fluxes are obtained from the entropy 
law; dfj^S^ > 0. Space-time evolution of the fluid is then obtained by solving energy- 
momentum conservation equations and relaxation equations for dissipative fluxes, 

a^T'''' = 0, (1) 
Dtt^"" = - — (tt'"' - 2j7V<''u">) - [u^n"^ + u''Tr^'^]Dux (2) 

Ttt 

DIl = - — {Il + Cd^u>') (3) 
Tn 



Eq.l is the conservation equation for the energy-momentum tensor, T'^" = {e+p)u^^u'^— 
[p + n)g'"^ + TT^'', e, p and u being the energy density, hydrostatic pressure and fluid 
4-velocity respectively. In Eq.2, D = uf^d^ is the convective time derivative, X/'^^u"^ = 
^{'V'^u'^ + 'Vu'^) - l{d.u){g'^'' - u^u") is a symmetric traceless tensor. Eq.2 and 3 are the 
relaxation equations for the shear stress tensor tt'"^ and bulk pressure H, with Tj^ and rjj 
as the relaxation times, rj and C are shear and bulk viscosity coefficients. We note that re- 
laxation equations for dissipative fluxes are non-trivial and may contain additional terms 
[16,17]. Eqs.1-3 are closed only with an equation of state (EOS) p = p{s). Any phase 
transition in the system can be incorporated through the EOS. Boost- invariant version 
of Eqs.l- 3 has been solved by many authors [1-9] ^ . Solution require initial conditions 
e.g. the initial or the thcrmalisation time, initial energy density, fluid velocity distribu- 
tion. In viscous hydrodynamics additionally, shear stress tensors and bulk pressure has 
to be initialised. A priori, they are unknown and only way to fix them is to confront the 
theory with experimental data. Fitting large number of parameters with experimental 
data is a complex process. To reduce free parameters, some of the parameters are fixed, 
e.g. initial transverse energy density/fluid velocity profile is fixed from Glauber model 
or color glass condensate model, initial shear stress tensors are fixed cither to zero or 
boost-invariant values etc. see [1-9]. All the parameters are not independent either. For 
example, initial time and central energy density are inversely related. Viscosity fiattens 
Pt spectra, an effect that can be mimicked by initial fluid velocity. Hydrodynamic models 
also require a freeze-out prescription. In ideal hydrodynamics, one generally use sudden 
freeze-out, where on some assumed freeze-out surface hydrodynamically evolving fluid 
suddenly transforms in to free-streaming particles. Most of the viscous hydrodynamics 
also uses the same algorithm. Recently, Dusling and Teaney [9] implemented dynam- 
ical freeze-out through the freeze-out condition, Treidu « 0.5. It satisfies the viscous 
hydrodynamics requirement that microscopic relaxation time is much smaller than the 
macroscopic inverse expansion rate. 



^ In [9] relaxation equations axe obtained in Ottinger and Gremala formulation [18]. 
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3. Hydro dynamical model's estimates oi rj/s from experimental data 

As mentioned earlier, one of the aims of causal hydrodynamics is to extract QGP 
viscosity from experimental data. In [8], charged particles multiplicity, radial flow and 
elliptic flow in i/SiVAr=200 GeV Au+Au collisions were analysed to obtain an estimate of 
rj/s. They also estimated systematic uncertainty both in theory and experiments, 77/s = 
0.1 ± 0.1{theory) ± 0.08{experiment). Considering several deficiencies of the model, it 
was concluded that QGP viscosity, r]/s > 0.5 can be excluded. A similar conclusion 
was drawn in [5] from the analysis of elliptic flow data. In the following, we discuss in 
brief some recent estimates of rj/s from <p meson data, centrality dependence and scaling 
violation of elliptic flow. 



Fig. 1. (color online) (left panel) Vari- 
ation of "^^^ and (Py) with central 
energy density in 0-5% Au-|-Au colli- 
sions, (right panel) centrality depen- 
dence of (p meson multiplicity, inte- 
grated V2 and mean px are compared 
with hydrodynamical simulations with 
central energy density as in table. 1. Ini- 
tial time Ti=0.6 fm, freeze-out temper- 
ature Tp=150 MeV. 




ri/s 





0.08 


0.12 


0.16 




35.5 ± 5.0 


29.1 ±3.6 


25.6 ±4.0 


20.8 ± 2.7 


n (MeV) 


377.0 ± 13.7 


359.1 ± 11.5 


348.0 ± 14.3 


330.5 ± 11.3 



Table 1 

Initial energy central energy 
density (e^) and temperature 
(Ti) from fit to (p multiplicity in 
0-5% Au+Au collisions. 



3.1. star's (j) meson data and QGP viscosity 

(f) mesons have several unique features, e.g. hidden strangeness, both hadronic and 
leptonic decay, not affected by resonance decays, mass and width are not modified in a 
medium etc. which make them an ideal probe to investigate medium properties in heavy 
ion collisions. We have analysed STAR data [21] on centrality dependence of (j) mesons 
multiplicity, mean pT and (integrated) V2 and obtained an accurate estimate of QGP 
viscosity over entropy ratio. Details can be found in [22]. With a lattice based EOS, 
where the confinement-deconfinement transition is a cross-over at Tc— 196 MeV, hydro- 
dynamic equations are solved with the code AZHYDRO-KOLKATA [1~3]. Variation of 
(j) multiplicity in 0-5% Au+Au collision, with the central energy density of fluid is shown 
in the left panel of Fig.l. For 7^/5=0-0.16, experimental multiplicity (the shaded region) 
can be fitted by changing the initial energy density. In tabid the central energy density 
and temperature required to reproduce the experimental multiplicity data are noted. 
More viscous fluid require less initial energy density /temperature as entropy is generated 
during the evolution. One note that with increasing viscosity mean pT increases and si- 
multaneous fit to experimental multiplicity and mean px is possible only for rj/s <0.12. 
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Fit obtained to the STAR data for centrality dependence of ^^3^, (Pt) ^^'^ ^2, from 
hydrodynamic evolution of fluid with viscosity to entropy ratio ry/s=0-0.16 are shown 
in the right panel, analysis indicate that the best fit to the data is obtained with 
ri/s = 0.07 ± 0.03 [22]. The central value is very close to the KSS bound, 77/s = l/An 
[23]. In [22], systematic uncertainty in 77/s due to uncertain initial conditions was also 
studied. Systematic uncertainty due to uncertain initial time Ti=0. 2-1.0 fm, freeze-out 
temperature Tf=140-160 MeV, percentage of hard scattering contribution in initial en- 
ergy density /=0-95%, initial transverse velocity Vr = tanh(ar), a=0-0.6, inaccuracy in 
the hydrodynamical code were evaluated to obtain, r]/s = 0.07±0.03(stot.) ±0.14(s?/s.). 
We note here that all the sources of systematic error e.g. uncertainty in initial stress 
tensors, relaxation time etc. are not included. Their inclusion will increase systematic 
uncertainty even more. 

3.2. Scaling property of elliptic flow and QGP viscosity 

Elliptic flow in non-central collisions is a key observable in establishing fluid like behav- 
ior of the medium produced in Au+Au collisions at RHIC. Scaling properties of elliptic 
flow (U2) has been studied by the STAR and PHENIX collaboration [24,25]. Both the 
collaborations observed 'universal scaling', W2/ ("■5(^2)0/1), elliptic flow scaled by the con- 
stituent quark numbers and charged particles V2 , of different particle species in different 
collision centrality, scales with KEx/nq, the transverse kinetic energy per constituent 
quark number. Initial eccentricity scaling or the constituent quark number scaling are 
not indicated in the data. 



In Fig. 2, we have studied universal 
scaling of elliptic flow for vr^, A'+, pro- 
ton, (j) and n in 0-60% Au-|-Au coUi- 
sions [26] . Only a few collision central- 
ities are shown for clarity. Contrary to 
experiments, in hydrodynamic simula- 
tions, irrespective of fluid viscosity, el- 
liptic flow does not follow the univer- 
sal scaling. The solid line in Fig. 2 is 
the approximate scaling function as ob- 
tained in the PHENIX experiment [25] 
(we have ignored the fluctuations of the 
scaling function). 

To measure the departure of simulated flows from the experimental 'universal scaling 
function', we defined a scaling violation function F, 

F = m^tcaled) - vTscaledf (4) 

In 0-60% Au-fAu collisions, simulated flows for tt^, A'+, p, (p and show minimum 
departure from universal scaling for rj/ s = 0.12 ± 0.03 [26]. 

In [26] initial eccentricity or constituent quark number scaling was also studied. Initial 
eccentricity scaling, or constituent quark number scaling is only approximate in ideal and 
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Fig. 2. (color online) Universal scaling in viscous hy- 
drodynamics. 
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viscous flow, though it appear that scaUng property is better satisfied in viscous than 
ideal fluid flow. A similar trend can be observed in Fig. 2. The result is contradictory to 
our expectation. Since viscous correction introduces a microscopic scale, it is expected 
that any scaling property that would have been observed in ideal fluid, will worsens 
in viscous fluid. The issue is discussed in [26]. In viscous evolution, elliptic flow has a 
contribution from the non-equilibrium part of the distribution function. In hydrodynamic 
evolution with lattice based EOS with cross-over transition, the viscous fluxes decreases 
rapidly. At freeze-out fluid behave more like an ideal fluid. Viscosity however changes the 
freeze-out surface. Changed freeze-out surface, but small viscous correction, may possibly 
be the reason for obtaining better scaling in viscous dynamics than in ideal dynamics. 

Table 2 



Best fitted r\l s as a function of collision centrality. 



coll. centrality 


0-10% 


10-20% 


20-30% 


30-40% 


40-50% 


50-60% 


0-60% 


rils 


± 0.03 


0.051 ± 0.008 


0.087 ± 0.004 


0.109 ± 0.003 


0.134 ±0.004 


0.169 ±0.005 


0.12 ± 0.01 



3.3. Centrality dependence of elliptic flow and QGP viscosity 

PHENIX measurements [27] for charged particles elliptic flow in 0-60% Au+Au colli- 
sions are shown in Fig. 3. In Fig. 3, the sohd, dashed, medium dashed and short dashed 
lines are AZHYDRO-KOLKATA simulations [28] for elliptic flow in fluid evolution with 
(i) rj/s—Q (ideal fluid), (ii) 77/5=6. 08, (iii) 77/5=0.12, and (iv) 77/5=0.16, respectively. In 
mid-central or in peripheral collisions, elliptic flow is over predicted in ideal fluid evo- 
lution. Data are better explained in viscous fluid evolution. It appear that elliptic flow 
data in more peripheral collisions demand more viscous fluid [28]. To be more quantative, 
from a analysis, we obtain the best fitted 77/5 in each collision centrality. They are 
listed in table. 2. In 0-60% collision centrality, -q/s varies between 0-0.17. 



Increase of viscosity with collision cen- 
trality is not contrary to the present 
paradigm that 77/s has a minimum, 
possibly with a cusp, around the crit- 
ical temperature T — Tc [29]. Rather 
it indicates the increasingly important 
role of hadronic matter in the develop- 
ment of elliptic flow in peripheral col- 
lisions. Viscosity to entropy ratio, as 
obtained from hydrodynamic analysis 
is averaged over the space-time. Both 
the QGP and the hadronic phase con- 
tribute to the average. Schematically, 
one can write. 



77 77 flQP 77 /^fl" 

j^{l~fHAD){-^) {Tqgp) + fHADi-^) (ThAd) (5) 




5 



where /had is the fraction of the hadronic matter, {^Y^^ {Tqgp) is the viscosity of QGP 
matter at average temperature Tqgp and {-)^°''^{Thad) is the viscosity of the hadronic 
matter at average temperature Thad- Detail analysis [28] indicate that compared to a 
central collision, in a peripheral collision, while Tqqp decreases, Thj^u remain approxi- 
mately same. Centrality dependence of extracted rj/s can be qualitatively understood as 
due to increased contribution of hadronic phase and decreased contribution of the QGP 
phase in peripheral collisions than in a central collisions. 
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Fig. 4. (color online) PHOBOS data for 
the centrality dependence of eccentric- 
ity scaled elliptic flow in Cu-|-Cu and 
Au-|-Au collisions at 62 and 200 

GeV. The solid line is the fit to the 
data with Eqs.6 and 8. We have used 
CT=3 mb, Ko=0.7 and = a/VS. The 
dashed line is the fit in the ideal fluid 
approximation. Quality of fit is com- 
paratively poor in ideal fluid approxi- 
mation. 
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4. Knudsen number, hydrodynamic limit for V2 and QGP viscosity 

Ideal hydrodynamics require complete thermalisation. The requirement is relaxed in 
dissipative hydrodynamics. Incomplete thermalisation in heavy ion collisions can be quan- 
tified in terms of Knudsen number, K = ^, Ris the characteristic size of the system and 
A is the mean free path. Applicability of hydrodynamics require that >> 1. In the 
Knudsen regime << 1 and hydrodynamics is inapplicable. The simple formula, 

proposed in [31] give qualitatively correct behavior of the experimental elliptic flow. In 
the limit of small Knudsen number, experimental flow approaches the ideal hydrodynamic 
limit — Uj'*, with a correction linear in K. In the other extreme limit of large 

oc K. In Eq.6, K^^ is a number of the order of unity. In [31], Knudsen number was 
expressed in terms of experimental multiplicity ^ , 

1 1 dN 

where S* is a measure of the transverse area of the collision zone, Cs is the square speed 
of sound of the medium and a is the inter-particle cross section. Eq.6 and 7 connect 
two experimental observables, elliptic flow and particle multiplicity and can be used to 
determine ideal hydrodynamic limit of elliptic flow (^) ' , the combination of parameters 
KqctCs- They have been used to estimate QGP viscosity from experimental data also. 
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The estimates vary between 7//s w 0.11 — 0.27 [32-34]. However, Eq.7 assume isentropic 
expansion and is invalidated in viscous evolution where entropy is generated. Recently in 
[37], Eq.7 was generalised to include the effect of entropy generation in viscous evolution. 



1 



1 dN 
S dy 



3Tm 



1 -3 



(8) 



One immediately observes that neglect of entropy generation during evolution will over 
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estimate K ^,hy the factor 



3TiTi 



(!) 



. Eqs.6 and 8 can be used to fit experimental 
data on particle multiplicity and V2- In Fig. 4, fit obtained to the PHOBOS data [35,36] 
with {ifY'' = 0.33 ± 0.12, (t;^! ) = 0.83 ± 0.51 is shown. Large uncertainty in (2^)'^ 

and (^7^7^ reflect the large uncertainty in PHOBOS data. Using simulation results 

of [22], ^7^!^ can be converted into more comprehensible viscosity to entropy ratio, 
7?/s = 0.17 ±0.10 ±0.20. 
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Fig. 5. (color online) For a parametric bulk viscosity (see left panel) AZHYDRO-KOLKATA simulations 
for pt spectra (middle panel) and V2 (right panel) in b=6.5 fm Au+Au collisions. r;/s=0.08. 



5. Bulk viscosity in heavy ion collisions 

Estimates of t]/s discussed above ignores the bulk viscosity. In general bulk viscosity is 
much smaller than shear viscosity. However, in QCD near the critical point bulk viscosity 
can be large [19,20]. In [6] effect of bulk viscosity on particle spectra and elliptic flow has 
been studied in detail. Strong growth of bulk viscosity near the critical point is offset 
by critical slowing down of the dynamics of bulk pressure, diminishing the bulk pressure 
contribution to spectra or flow [6] . Recently, we have included the effect of bulk viscosity 
in AZHYDRO-KOLKATA. Detail resuhs will be published later. In Fig.5, AZHYDRO- 
KOLKATA simulations for px spectra and V2 in b=6.5 fm Au+Au collisions are shown. 
Pt spectra and V2 are marginally changed when effects of bulk viscosity is included. 
Indeed, the sensitive observable V2 also show little modification. It is understood. Bulk 
viscosity do not introduce additional asymmetry in the system and elliptic flow is largely 
unaffected by bulk viscosity. Estimates oi rj/s obtained by neglecting bulk viscosity will 
not be altered significantly with inclusion of bulk viscosity. 
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6. Summary 

To summarise, we have briefly discussed causal dissipative hydrodynamics. Some recent 
attempts to extract QGP viscosity from comparison of hydrodynamical simulations with 
experimental data arc also discussed. Extracted values of viscosity to entropy ratio, from 
a variety of experimental data, lies within a narrow range, rj/ s «0. 08-0. 20. However, due 
to uncertain initial conditions in hydrodynamic simulations, systematic error in extracted 
values of -q/s could be large, ^200% or more. It appear that precise estimate of QGP 
viscosity from experimental data will not be possible, unless inputs of hydrodynamical 
models are constrained. 
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